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SECTION  I  INTRODUCTION 

Structural  damage  to  submerged  shell  structures  excited  by  a  shock  wave  is 
an  inherently  nonlinear  phenomenon.  Fundamental  to  the  study  of  such  phenomena 
is  the  formulation  of  a  consistent  and  reliable  theory  to  predict  the  nonlinear 
transient  response  of  a  shell  structure.  Naturally,  a  complete  investigation 
of  the  problem  must  include  the  fluid-structure  iteration.  Here  only  the 
nonlinear  structural  response  will  be  considered. 

The  modal  superposition  method  has  for  many  years  been  the  dominant 
dynamic  analysis  method  for  linear  systems .  Recently,  it  has  been  suggested 
[Ref.  1,2,3]  that  the  modal  superposition  method  can  be  successfully  applied  to 
nonlinear  problems,  it  would  appear  that  modal  superposition  is  competitive 
with  direct  integration  for  many  nonlinear  problems  and  may  be  better  for 
problems  involving  "wide-band"  governing  equations. 

It  has  also  been  recognized  that  the  modal  superposition  method  is  aided 
when  it  is  supplemented  with  a  static  solution  calculated  for  the  load  applied 
at  particular  times.  In  dynamics  treatises  of  20  or  30  years  ago,  this  method 
was  referred  to  by  a  variety  of  names,  e.g.,  the  acceleration  method.  It 
is  now  sometimes  called  the  residue  flexibility  method.  This  method  was 
recognized  as  improving  convergence,  at  least  during  the  period  for  which  the 
loading  is  applied,  it  is  now  believed  [Ref.  4]  that  the  use  of  the  static 
response  will  account  for  the  effect  of  higher  modes. 
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SECTION  II  FORMULATION  OF  THE  METHOD 


1.  DEVELOPMENT  OF  THE  METHOD 


It  is  proposed  that  a  modal  superposition  analysis  be  employed 
for  nonlinear  problems.  The  static  response  at  each  point  in  time 
is  to  be  used  to  supplement  the  modal  susnation.  The  nonlinearities 
can  be  incorporated  using  an  incremental  or  a  pseudo-force  approach. 

For  the  formulatioin  to  be  outlined  here,  use  of  the  pseudo-force  method 
is  probably  the  most  straight  forward  technique. 

Let  the  p  normal  mode  shapes  for  the  structure  at  time  t  ■  o  be 


{$  > 


{ 4  } 
2 


-  t*  > 
P 


(1) 


These  should  be  computed  at  the  outset  of  the  solution. 

At  each  time  t  calculate  the  static  displacements  due  to  the 
load  occurring  at  time  t.  use  the  updated  stiffness  characteristics. 
Let  (u  )  be  this  solution,  i.e.  the  solution  to  the  static  problem 

_ x.  x _ 


Let  {u}  be  the  physical  degrees  of  freedom.  Then,  set 

(u>  «  U  ‘  u  1  <q>  -  CR]{q)  (2) 

NLi 


Where 


m  -  t  C4>  (<j> 2>  -  .  .  ) 

(3) 

(q)  “  modal  coordinates 

[R]  -  C  ♦  u  ] 

NL 

(4) 

The  modal  equations  would  then  be  coupled 

[R)T(m] [R] {q>  +  [RMkHRHq)  -  [RjV)  (5) 
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where 


[m]  is  the  mass  matrix 
[k]  is  the  stiffness  matrix 
{F}  is  the  loading  vector 

Although  Eqs.  (S)  are  coupled,  note  that 

[R]T[m][R]  and  [R]Tfk][R] 

are  of  much  smaller  order  than  [m]  and  [k],  In  fact,  for  linear  systems  they 
take  the  special  form 


For  some  cases,  an  efficient  algorithm  for  solving 
Eq.  (5)  can  probably  be  established. 


As  an  alternative  to  this  procedure,  the  {u_ )  can  be 
transformed  to  a  more  favorable  set  of  coordinates.  Compute  {u^} 
from  (u^)  which  is  orthogonal  to  .  .  .  ,  p)  •  Set 

{u}  -  *.  i  ](q)  -  [R]  lq }  (6) 

NL 


The  modal  equations  will  be  uncoupled  here.  However,  in  order  to 
achieve  the  uncoupled  modal  equations  it  is  necessary  to  calculate 
from  tu^,. 


2.  DISCUSSION 


The  proposed  procedure  is  a  nonlinear  version  of  the  residual 
flexibility  or  acceleration  approach  for  modal  superposition.  The 
advantage  in  this  technique  over  the  use  of  updated  { <f> .  }  may  lie  xn 
both  greater  efficiency  and  accuracy.  There  is  some  question  as  to 
the  effectiveness  of  the  use  of  updated  mode  shapes,  say  0^},  since 
{?}  will  be  formed  of  linear  combinations  of 

To  achieve  the  desired  efficiency  and  accuracy,  it  is  proposed 
that  a  novel  static  reanalysis  procedure  be  employed  to  update  the 
static  solutions  at  time  t.  This  is  the  key  to  the  solution 
technique . 


3.  STATIC  REANALYSIS  METHODS 


The  success  of  the  proposed  method  depends  on  the  availability  of 
an  efficient  and  reliable  method  for  the  static  reanalysis  of  a  structural 
system.  Before  developing  such  a  method  we  will  review  some  of  the 
currently  available  reanalysis  techniques.  The  goal  of  reanalysis 
technology  is  to  analyze  a  modified  structure  by  taking  advantage  of  in¬ 
formation  from  the  solution  to  the  original  or  unmodified  structure  to 
avoid  a  complete  analysis  of  the  modified  structure. 

The  methods  will  be  grouped  as  approximate  and  exact  techniques. 

The  approximate  methods  can  be  categorized  as:  perturbation  (Taylor  series) 
method,  iteration  method,  and  reduced  basis  method.  These  are  generally 
applicable  to  small  magnitude  modifications  distributed  over  a  large 
portion  of  the  structural  system.  In  constrast,  the  exact  methods  are 
applicable  to  localized  modifications  with  arbitrary  magnitudes.  Application 
of  the  exact  method  requires  that  the  structural  configuration  and  the 
validity  of  the  finite  element  model  are  not  affected  by  the  modification. 

The  underlying  principles  of  the  exact  methods  involve  the  concep¬ 
tion  of  parallel  elements,  initial  stress,  and  pseudo  loads.  Pure 
algebraic  formulas  such  as  the  Sheraan-Morrison  formula  (5)  and  its 
generalization,  the  Woodbury  formula  (5)  can  also  be  directly  used  to 
develop  exact  modification  methods.  In  this  section,  a  direct  method 
based  on  the  Woodbury  formula  due  to  Argyris  and  Roy  (6) ,  a  successive 
modification  method  proposed  by  Kavlie  and  Powell  (7)  and  the  pseudo 
load  method  developed  by  the  authors  are  discussed.  The  last  method 
although  similar  to  that  of  Argyris  and  Roy  is  derived  by  only  using  the 
system's  linearity. 

In  addition  to  reanalysis  the  pseudo  load  approach  provides  closed 
for  semi-closed)  form  relationships  between  static  response  variables 
and  the  design  parameters.  This  enables  the  design  (i.e.,  select  design 
parameters)  for  prescribed  responses. 


Notation 

Static  Reanalysis 

A* 

b  =  a  Boolean  matrix  which  expands  A K  into  AK, 

_  M 

AK  =  bT  AK  b 
B  =  system  bandwidth 
dt  *  ith  'esign  par /  ?eter 
D  =  number  '  J  d«_.  ign  variables 


r  **iWm*.  -  ■•  '^V^Vjw-  • 


K  =  system  stiffness  matrix 

K.  =  modified  system  stiffness  matrix  with  the  1  through 
ith  non-zero  columns  of  AK  included 

AK  =  stiffness  modification  matrix  of  system 


AK.  =  X.  = 
— 1  -l 


AK  with  all  null  columns  removed 


ith  non-null  column  of  AK 


AK  =  AK  with  all  null  rows  removed 


lower  triangular  matrix 


n  = 
c 


=  number  of  non-null  columns  (rows)  of  AK 


external  load  vector 


=  number  of  reduced  basis  vectors 
=  transformation  matrix 
=  original  system  solution 


modified  system  solution 


t(i)  , 


=  ith  basis  vector  used  in  the  reduced  basis  method 

=  modified  system  solution  after  the  1  through  ith 
non-null  columns  of  AK  are  included 


W<k>  = 

-l 


=  h'ii 

AK^ ,  ith  column  of  X 
ith  row  of 


Approximate  Methods 


Taylor  Series  Expansion  ( 8  ) 

The  static  responses  can  be  expanded  in  Taylor  series  as  functions 
of  the  design  parameters,  for  example, 

D  3  U  D  D  32  U 

U  =  U  +  I  -  Ad.  +  I  I  -  Ad.  Ad.  ns 

1  1  j  \>) 

i=l  3d.  i=l  j=l  3d. 3d. 

1  1  J 

+  higher  order  terms 

The  derivative  terms  are  evaluated  by  differentiating  the  system's 
static  equilibrium  equation, 

K  U  =  P  (8) 


assuming  that  the  external  load  vector's,  P,  change  is  negligible  when  a 
design  parameter  changes.  The  first  two  derivatives  are  expressed  by, 

3  U  .  3  K 

=  -  K  U  (9) 


32  U 

3d ! 3d . 
i  J 


32  K  3  K  3  U 

-  (3d.3d.  y+3d?  3d7 
1  J  1  J 


3  K  3  U 

+  3d~  3d~  * 
J  i 


(10) 


Note  that  in  practice  K  is  not  inverted,  instead  the  triangularized  form 
of  K  obtained  from  the  unmodified  system's  solution  is  used  to  solve  for 
the  vectors  in  equations  (  9  )  and  (10).  In  addition,  to  maximize  the 
computational  efficiency,  advantage  should  be  taken  of  the  sparsity  of 
the  derivatives  of  K.  Ordinarily  only  the  first  order  terms  of  equation 
(  7  )  are  retained. 

Storaasli  and  Sobieszczanski  (_9_)  have  applied  the  first  order 
Taylor  series  method  to  a  5000  dof  fuselage  midsection  model.  Less  than 
16%  error  was  observed  for  single  element  variation  of  -100%  (i.e. 
element  removal)  to  +500%,  and  for  -50%  to  50%  change  for  multi¬ 
elements.  The  incorporation  of  second  order  terms  would  undoubtedly 
increase  this  solution  accuracy,  but  at  the  expense  of  increased  com¬ 
putational  burden. 

Iteration  Method 

The  equilibrium  equation  of  the  modified  system  may  be  written  as, 

K  0  =  P-AK  0  (11) 

From  this  form  an  iterative  solution  approach  is  derived  in  the  follow¬ 
ing  recursive  form, 

K  U(l+l)  =  P  -  AK  U(l) 


(12) 


Back  substitution  is  employed  to  solve  for  U  1  ,  utilizing  the  tri- 

angularized  form  of  K  from  the  unmodified  system's  solution.  The  itera¬ 
tion  is  started  using, 

U(o)  =  U  =  unmodified  system  response 
Noor  and  Lowder  ( ]_)  suggested  a  method  comprised  of  the  Taylor  series 
and  iteration  formulas.  First  0^°^  is  computed  via  a  Taylor  series 


method,  i.e. 


D  9  U 


U  =  U  +  2 


I  Ad. 

.  ,  9d .  l 

i=l  i 


Then  IT  is  used  to  initiate  the  iterative  process  in  eq.(12),for  i>0. 
Usually  one  or  two  iteration  cycles  are  reported  sufficient  to  obtain 
good  results. 

Reduced  Basis  Method 

Fox  and  Miura  (10)  introduced  a  procedure  analogous  to  assumed 
modes  to  reduce  the  dimension  of  the  modified  structure's  static  equi¬ 
librium  problem.  A  set  of  known  system  deflection  vectors 
(r) 

...,  U  are  assumed  to  span  the  solution  space.  Since  r«n  =  number 
of  dof,  this  is  obviously  an  approximation.  The  approximate  solution  of 
the  modified  system  then  follows  as, 


U  =  l  C.  U 

"  j=l  J" 


K  C  =  P 


K  =  R  (K  +  AK)  R 
T 

P  =  R  P 


C  =  (Cj  C2 


(rxl) ,  R  = 


(rxl)  (16) 

(rxr)  (17) 

!  U(r)]  (nxr)  (18) 


Fox  and  Miura  indicated  this  is  equivalent  to  applying  the  Ritz-Galerkin 
principle  in  the  subspace  spanned  by  the  set  of  vectors  ...  U^r^. 

Although  suggestions  were  provided,  a  systematic  approach  to  calculating 
the  basis  vectors  and  deciding  the  value  of  r  was  not  indicated.  Noor 
and  Lowder  (_8J  suggested  the  use  of  the  following  as  basis  vectors, 

(i)  (i+D  a  u 

U  =  U  ,  U  =  jr-r—  i=l  ,2, .  .  .D  (19) 
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where  D  is  the  number  of  design  parameters.  Excellent  results  were 
obtained. 


Exact  Methods 

Successive  Modification  Methods  (SMM) 

This  method,  developed  by  Kavlie  and  Powell  (7),  incorporates  the 
modification  matrix,  AK,  a  single  non-null  column  at"  a  time.  The  basis 
of  this  approach  is  the  Sherman-Morrison  formula  (see  Householder  (5  )) 
for  the  inverse  of  a  modified  matrix,  i.e. 


(A  +  X  YT)_1  =  A"1  -  (A-1  X)(YT  A~ 1 ) / ( 1  +  YT  A"1  X)  (20) 

where  X  and  Y  are  column  vectors  representing  the  modification  to  A. 

The  system's  static  equilibrium  equations,  before  and  after  struc¬ 
tural  modifications  are, 

K  U  =  P,  (K  +  AK)  U  =  P  (nxl)  (21) 

Consider  the  special  case  when  just  column  j  of  M  is  non-null.  Desig- 

<v 

nate  this  column  as  X  =  AK  .  Then 

AK  =  X  YT  =  AK  YT,  YT  =  [0  ...  0  1  0  ...  0]  (22) 

column  j 

and  equation  (  20  )  implies 

(K  +  AK)"1  =  K_1  -  K'1  X  YT  K_1/(l  +  YT  K_1  X)  (23) 

To  apply  equation  (23)  directly  would  explicitly  require  K  .  Alter¬ 
natively,  post  multiplying  this  equation  by  P  yields 


U  =  U  -  W  YT  U/(l  +  YT  W),  K  W  =  X  (24) 

where  W  is  calculated  by  back  substitution  using  the  triangularized  form 
of  K  from  the  unmodified  system's  solution. 

For  the  more  general  case  place  the  n  non-null  columns  of  AK  in 
the  columns  of  the  nxn  matrix  X.  The  non-null  column  numbers  of  AK  are 
designated  by  the  set  of  ascending  integers, 


jl  j2 


(25) 


T 

With  these  definitions  the  ith  row  of  the  n  xn,  Boolean  matrix  Y  which 
expands  X  into  AK  is  null  except  for  the  number  one  (1)  in  column  j.. 
The  stiffness  modification  matrix  is  then  given  by, 


-v^f  •*  -.t-r-  >  •  -h-  &*&***!**&& 


AK  =  X  Y 

(nxn)  (nxn  )  (n  xn) 
c  c 


=  AK  Y1 


For  example  consider  a  10  dof  system  with  modifications  to  K  in  columns 
2,  5,  7  and  9,  then 

nc  “  Jj  =  2,  J2  -  5,  -  7,  -  9 

AK  =  [0  !  AK„  !  0  !  0  !  AKC  !  0  !  AK,  !  0  '  AK-  !  0] 

—  -  1  — 2  1  -  1  -  1  — 5  1  -  1  — 7  1  -  t  — 9  1  - 


X  =  [AK.  .  AK„  .  AK,  1  AK- ]  =  [X.  •  X  <  X-  1  X. ]  =  AK 

-  — 2  1  — 5  1  — 7  1  — 9  -1  1  -2  1  -3  1  -4'  — 


0100000000 


Ya  =  0000100000 


0000001000 


0000000010 


Then  if  more  than  one  column  of  K  is  modified,  the  successive  modifica¬ 
tion  method  of  Kavlie  and  Powell  is  employed.  This  assumes  the  follow¬ 
ing  recursive  form; 


U.  =  U.  -  W^1_1)  YT  U.  ./( 1  +  yt  wfI_1)) 

-x  -l-l  -1  -1  -l-l  -1-1 


W^k)  =  K71  X. 


U  =  U 


where  U.  is  the  system's  displacement  vector. after  the  first  i  non-null 
columns  of  AK  are  considered.  Note  that  W.  '  is  the  system’s  displace¬ 
ment  vector  at  the  stage  when  k  non-null 1  columns  of  AK  are  included, 
with  X.  as  the  external  load  vector.  These  vectors  are  calculated  by 
the  formulas 

w(k)  _  w(k-l)  _  w(k_1)  yj  W^'^/Cl  +  y!  wfk”1^)  (29> 

— 1  —1  -k  —  k  —1  — k  — k 


K  W.  =  X. 
- 1  -1 


i  =  1,2,  ...  n 


where  W.°  is  calculated  by  back  substitution  using  the  triangularized 
form  of1  K  from  the  unmodified  system's  solution.  Equations  (27) 
through  (30)  follow  directly  from  the  version  of  the  Sherman-Morrison 
formula  given  by  equation  (  24  ) .  A  flow  diagram  for  the  10  dof-n  =4 
system  is  shown  in  Figure  1. 


Kavlie  and  Powell  provide  operation  counts  for  this  method,  and 
emphasize  that  any  reanalysis  scheme  should  be  carefully  scrutinized  as 
to  when  it  is  more  efficient  than  complete  reanalysis. 

Direct  Modification  Method 

This  method  was  developed  by  Argyris,  et.  al.  (_6_) .  The  derivation 
.s  based  on  the  Woodbury  formula  (ref.,  5 ) , 

(A  +  E  C  FV1  =  A-1  -  A"1  E  (O'1  +  FT  A-1  E)'1  FT  A-1  (31) 


however  a  derivation  from  physical  arguments  is  also  included.  For 
reanalysis  of  the  static  equilibrium  equation  substitute 


A  =  K,  E  =  bT,  F  =  bT,  C  =  M 

(32) 

AK=bTAKb=ECFT 

(33) 

T 

into  equation  (  31  )  where  b  =  Y  (see  equation  (  26  )). 

This  yields 

(K  +  AK)-1  =  (K  +  bT  AK  b)'1  =  K_1  -  K_1  bT  (AK-1  +  b  K_ 

1  bV'b  K'1 
(34) 

Post  multiply  equation  (34)  by  P  and  use  equation  (  21 ) 

.  This  gives 

•V/ 

AU  =  0  -  U  =  -  (LT)_1  Z  2'1  (g"1  +  AK)'1  AK  b  U 

(35) 

where  from  the  unmodified  system  solution  the  Cholesky  decomposition 
of  K  is, 

T 

K  =  L  Ll 

(36) 

and  L  Z  =  bT  ,  §  =  ZT  Z 

(37) 

Computationally,  the  following  procedure  was  suggested, 


(1)  Compute  r’=  M  b  U,  (2)  Evaluate  Q  and  Z, 

T  -1 

(3)  Triangularize  Q;  Q  =  Lq  ^  and  evaluate  Q  , 

(4)  Triangularize  (Q  1  +  AK) ;  (Q  1  +  AK  =  L  1-3, 

s  s 

~  T  T  ~ 

(5)  Compute  g  and  r”  from  Lg  Ls  g  =  r’and  Lq  Lq  r’ ’  =  g, 

(6)  Compute  r ”  ’  =  Z  r”,  finally 

(7)  Solve  LT  (0  -  U)  =  -  r”’. 

Argyris  emphasized  that  the  preceeding  algorithm  is  more  efficient 
than  complete  reanalysis  only  when  n  is  approximately  less  than  0.75B. 
Also  mentioned  is  that  Argyris' s  approach  is  generally  more  efficient 
than  Kavlie  and  Powell’s  as  expressed  by  equation  (  27  )- (  30  ) . 
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A  New  Method  Pseudo  Load  Method 


It  is  hypothesized  here  that  the  modifications  actually  behave  as 
"closed- loop"  deflection  dependent  loadings,  forcing  the  unmodified  system. 
Furthermore,  suppose  that  for  a  localized  modification,  the  added  pseudo  load¬ 
ing  only  is  dependent  on  *  small  portion  of  the  displacements.  The  modified 
system  solution  employs  the  available  information  from  the  solution  of  the 
unmodified  system.  This  permits  the  modified  system  dispalcements  to  be 
expressed  as  functions  of  the  displacements  associated  with  only  that  portion 
of  the  system  where  modifications  occur.  Partitioning  the  modified  system 
displacement  vector  and  equating  the  modification  related  displacements  yields 
a  reduced  set  of  equations.  The  following  analytically  summarizes  the  pro¬ 
ceeding  concepts. 


Substitute  equation  (  26  )  into  equation  (  11  ) 


KU  =  P-  XYTU 
U  =  U  -  K-1  X  YT  U 


Employ  the  previous  notation  of  equation  (  24  ) . 
expressed  as 


0  =  u  -  w  yt  u 


T 

The  Boolean  matrix  Y  condenses  U  to  the  form, 


-  T  - 
U  =  Y  U 


(nxl)  (38) 

(nxl)  (39) 

Equation  ( 39  )  is 


(nxl) 


(ncxl) 


This  condensed  displacement  vector  containsonly  the  coordinates  which 
correspond  to  non-null  columns  of  AK,  i.e.  U.  ,  i=l,2,...n  .  Equation 
(  40 )  becomes ,  ^i 


U  =  U  -  W  U 


(nxl) 


By  shifting  rows  ( j . ,  i=l,2,...,n  )  in  U,  U  and  W  to  the  top  in 
equation  (42)*  it  is  seen  that 


n  xn 
c  c 


(nxl) 


(n-nc) 


(n-nc)  (n-nc) 


The  condensed  displacement  vector  is  then  calculable  from 


(I  +  W)  U  =  U 
c 


(ncxD 


(44) 


After  which  the  remaining  displacements  are  found  using 


U  =  U  -  W  U 

-r  -r  -r  - 


(n-nc)xl 


(45) 


The  computation  of  W  in  equation  (  40  )  proceeds  by  solving  the  separate 
sets  of  equations  K  W.  =  X. ,  i=l...n  .  These  sets  are  solved  by  back 
substitution  utilizing1  the  Vriangularized  form  of  K  from  the  unmodified 
system's  solution. 

The  pseudo  load  method  appears  to  resemble  the  direct  modification 
method.  Compared  with  the  successive  modification  method,  this  method 
treats  all  modifications  simultaneously.  In  addition,  the  derivation  of 
the  pseudo  load  method  is  simple  and  direct.  The  unmodified  system's 
displacement  vector  and  triangularized  form  of  K  may  be  obtained  from 
any  appropriate  finite  element  code.  A  simple  truss  problem  is  included 
to  illustrate  the  pseudo  load  method. 

The  geometric  and  material  properties  of  the  truss  members  are 
shown  in  Figure  2  .  Employing  the  direct  stiffness  method  (ref.  JL1)  with 
the  appropriate  boundary  conditions  provides  the  static  equilibrium 
equation, 


5  10  0 

V 

0 

0.125  x  106 

15  0-4 

U2 

= 

1000 

0  0  5  -1 

U3 

0 

f 

o 

1 

1 

*— » 

-V 

0 

Cholesky  decomposition  of  K 

is 

K  = 

T 

L  L  ; 

K  U  =  P 


(46) 
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Figure  2.  Truss  used  to  illustrate  the  pseudo  load  method  of 
static  reanalysis 
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The  solution  for  U  proceeds  in  the  standard  fashion, 

P  =  LLTUi=>P  =  LY  ,  y  =  ltu 

U  =  (2/11)  x  10'3  [-6  30  5  25]  (*7> 

Now  suppose  that  the  cross-sectional  area  of  member  4  is  changed  so  that 
A4  «=>A’  =  /2/2  (1  +  a).  The  responses  of  the  modified  system  due  to 

P,  as  a  function  of  a,  can  be  calculated  by  the  pseudo  load  method.  The 
required  modification  matrices  and  indices  for  the  problem  are, 


0 

3  0  0  “ 

0  0 

0  0 

0 

3  0  0 

, 

X  = 

0  0 

,  Y= 

0  0 

AK  = 

0 

3  e  -e 

e  -e 

1  0 

0 

3  -e  e 

-e  e 

_°  :_ 

u„ 

U, 

3 

1 

0  = 

_C4. 

,  o  = 

-r 

1 

CM 

o 

_ 1 

e  =  0. 

125 

<  106  •  a 

»  n 

CM 

It 

U 

j  1  -  3 

j2  = 

Using  the  decomposed  form  of  K,  the  matrix  W  is, 


-1 

5 

-1 


I  W  -  [  W 


-1 


W2]  =  - 
*  11 


1« 


-6 


6 


The  modified  system  responses  are  from  equations  ( 44  )  and  (45  ) 


U3  =  ((1  +  6oc/il)U3  +  a  U4/ll)/(l  +  7«/ll)  (48> 

U4  =  (6<xU3/11  +  (1  +  °c/ll)U4)/(l  +  70C/11)  (49) 

Uj  =  Uj-«(U3  -  U4)/U  (50) 

U2  =  u2  "  5a(-°3  +  V/U  (51) 

Results  like  equations  (48  )  and  (  49  )  are  generally  very  useful  since 
they  provide  (1)  an  efficient  basis  for  performing  a  parametric  study  of 
the  static  response,  and  (2)  means  for  prescribed  static  ^responses.  The 
second  use  is  illustrated  by  prescribing  either  U„  or  U4  in  equations 
(48)  or  (49  )  and  then  calculating  the  required  iactor  a,  i.e.  from 
equation  (  48  )  , 

«  =  11 (U3  -  U3)/(7  U3  -  6  U3  -  U4)  (52) 

For  example  if  it  is  desired  for  U_  =  0.9  !L,  then  from  equation  (52), 
a  =  -0.234,  and  from  equations  (  48;- (51  ), 

U  =  [-  6.5  32.5  4.5  28.0]  x  (2/11)  x  10~3 

If  from  equation  (52),  «  <  -  1.0,  then  the  prescribed  value  of  U~  is 
not  realizable  by  only  modifying  member  4.  i 
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4.  CONCLUSIONS  AND  RECOMMENDATIONS 

A  new  method  for  the  nonlinear  transient  response  of  structures,  in 
particular,  submerged  shell  structures,  is  formulated.  This  is  a  modal 
superposition  technique  using  residue  flexibility  in  which  the  static 
response  at  each  time  is  computed  using  an  efficient  static  reanalysis 
technique. 

This  preliminary  development  of  the  method  appears  to  indicate  that 
this  approach  holds  promise.  Some  numerical  tests,  which  are  not  reported 
here,  confirm  this.  The  method  seems  to  deserve  further  development. 
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DEPARTMENT  OF  ENERGY  CONTRACTORS 


Lawrence  Livermore  National  Lab 
ATTN:  S.  Erickson 


Deputy  Chief  of  Staff 
Research,  Development,  A  Acq 
Department  of  the  Air  Force 
ATTN:  R.  Steere 
ATTN:  AFRDQI 

Deputy  Chief  of  Staff 
Logistics  A  Engineering 
Department  of  the  Air  Force 
ATTN:  LEEE 

Foreign  Technology  Div 
Air  Force  Systems  Command 


ATTN: 

SD6G 

ATTN: 

S06F,  S.  Spring 

ATTN: 

TQTD 

ATTN: 

NIIS,  Lib 

United  States  Air  Force  Academy 
ATTN:  DFCEM,  W.  Fluhr 


Los  Alamos  National  Lab 

ATTN:  M/S634 ,  T.  Dowler 
ATTN:  MS  670,  J.  Hopkins 
ATTN:  R.  Whitaker 
ATTN:  R.  Sanford 
ATTN:  A.  Davis 
ATTN:  Reports  Lib 
ATTN:  G.  Spillman 

Oak  Ridge  National  Laboratory 
ATTN:  Civil  Def  Res  Proj 
ATTN:  Central  Rsch  Lib 

Sandia  National  Lab 
Livermore  Lab 

ATTN:  Library  A  Security  Classification  Div 

Sandia  National  Lab 
ATTN:  3141 
ATTN :  L .  Vortman 
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DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 


BDM  Corp 

ATTN:  T.  Neighbors 
ATTN:  A.  Lavagnino 
ATTN:  Corporate  Lib 

Bel  1  Tel ephone  Labs 

ATTN:  J.  White 

California  Institute  of  Technology 
ATTN:  T.  Ahrens 


NKF  Engineering  Associates,  Inc 
ATTN:  R.  Belsheim 

Pacific-Sierra  Research  Corp 
ATTN:  H.  Brode 

Pacifica  Technology 
ATTN:  R.  Bjork 
ATTN:  G.  Kent 


Physics  Applications,  Inc 

Cal i fornia 

Research  A  Tchnology,  Inc 

ATTN 

C.  Vincent 

ATTN 

Lib 

ATTN 

K.  Kreyenhagen 

Physics  International  Co 

ATTN 

S.  Schuster 

ATTN 

J.  Thomsen 

ATTN 

M.  Rosenblatt 

ATTN 

F.  Sauer 

ATTN 

L.  Behrmann 

Columbia  University 

ATTN 

E.  Moore 

ATTN 

F.  DiMaggio 

ATTN 

Tech  Lib 

ATTN 

H.  Bleich 

University  of  Pittsburgh 

University  of  Denver 

ATTN 

M.  Williams,  Jr 

ATTN 

Sec  Officer  for  J.  Wisotski 

Science  Applications,  Inc 

Electric  Power  Research  Institute 

ATTN 

Tech  Lib 

ATTN 

G.  SI  iter 

Southwest  Research  Institute 

General  Dynamics  Corp 

ATTN 

W.  Baker 

ATTN 

0.  Mador 

ATTN 

A.  Wenzel 

ATTN 

M.  Pakstys 

ATTN 

J.  Miller 

SRI  International 

ATTN 

G.  Abrahamson 

Kaman  AviDyne 

ATTN 

W.  Wilkinson 

ATTN 

N.  Hobbs 

ATTN 

A.  Florence 

ATTN 

Lib 

ATTN 

G.  Zartarian 

Systems,  Science  and  Software 

ATTN 

R.  Ruetenik 

ATTN 

K.  Pyatt 

ATTN 

R.  Sedgewick 

Kaman  Sciences  Corp 

ATTN 

Lib 

ATTN 

F.  Shelton 

ATTN 

D.  Grine 

ATTN 

Lib 

ATTN 

T.  Riney 

ATTN 

T.  Cherry 

Kaman  Sciences  Corp 

ATTN 

T.  McKinley 

ATTN 

D.  Sachs 

Teledyne  Brown  Engineering 

Kaman  Tempo 

ATTN 

J.  Ravenscraft 

ATTN 

DASIAC 

Tetra  Tech 

Inc 

Karagozian 

and  Case 

ATTN 

L.  Hwang 

ATTN 

J.  Karagozian 

TRW  Defense  &  Space  Sys  Group 

Lockheed  Missiles  A  Space  Co,  Inc 

ATTN 

Tech  Info  Ctr 

ATTN 

T.  Geers 

ATTN 

W.  Llpner 

ATTN 

B.  Almroth 

ATTN 

D.  Jortner 

ATTN 

Tech  Info  Ctr 

ATTN 

B.  Sussholtz 

ATTN 

A.  Narevsky 

Lockheed  Missiles  A  Space  Co,  Inc 

ATTN 

A.  Feldnan 

ATTN 

TIC-Library 

ATTN 

P.  Bhutta 

MAT  Company 

TRW  Defense  &  Space  Sys  Group 

ATTN 

D.  McNaight 

ATTN 

P.  Dai 

ATTN 

F.  Pieper 

Management  Science  Associates 

ATTN 

E.  Wong 

ATTN 

K.  Kaplan 

ATTN 

G.  Hulcher 

McDonnell  Douglas  Corp  Weldl lager  Assoc,  Consulting  Engrg 

ATTN:  R.  Halprin  ATTN:  J.  Isenberg 


26 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 


University  of  Virginia 
4  cy  ATTN:  T.  Giacofci 
4  cy  ATTN:  W.  Pi  1  key 


Weidlinger  Assoc,  Consulting  Engrg 
ATTN:  J.  McCormick 
ATTN:  M.  Baron 


